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We discuss the four-point correlation function, or the trispectrum in Fourier space, of CMB 
temperature and polarization anisotropics due to the weak gravitational lensing effect by interven- 
ing large scale structure. We discuss the squared temperature power spectrum as a probe of this 
trispectrum and, more importantly, as an observational approach to extracting the power spectrum 
of the deflection angle associated with the weak gravitational lensing effect on the CMB. We ex- 
tend previous discussions on the trispectrum and associated weak lensing reconstruction from CMB 
data by calculating non-Gaussian noise contributions, beyond the previously discussed dominant 
Gaussian noise. Non-Gaussian noise contributions are generated by lensing itself and by the cor- 
relation between the lensing effect and other foreground secondary anisotropies in the CMB such 
■ as the Sunyaev-Zel'dovich (SZ) effect. When the SZ effect is removed from temperature maps us- 

f*"^ | ing its spectral dependence, we find these additional non-Gaussian noise contributions to be an 

order of magnitude lower than the dominant Gaussian noise. If the noise-bias due to the domi- 
nant Gaussian part of the temperature squared power spectrum is removed, then these additional 
J_( , non-Gaussian contributions provide the limiting noise level for the lensing reconstruction. The tem- 

Oh perature squared power spectrum allows a high signal-to-noise extraction of the lensing deflections 

' and a confusion-free separation of the curl (or B-mode) polarization due to inflationary gravitational 

waves from that due to lensed gradient (or E-mode) polarization. The small angular scale temper- 
ature and polarization anisotropy measurements provide a novel approach to weak lensing studies, 
I | complementing the approach based on galaxy ellipticities. 

> : 
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\Q . I. INTRODUCTION 

O ■ 

The well understood features in the angular power spectrum of cosmic microwave background (CMB) anisotropies, 
such as acoustic peaks and the damping tail jlj, allow a useful probe of the cosmology J^|. The ability of CMB 
anisotropies to constrain most, or certain combinations of, parameters that define cold dark matter cosmological 
models with adiabatic initial conditions has driven a significant number of experiments from ground and space. These 
experiments include NASA's MAP mission^] and, in the long term, ESA's Planck surveyor^. With the successful 
detection of the first two or three acoustic peaks at degree angular scales [|| , the anisotropy experiments have started 
to focus on fluctuations on angular scales of arcminutes or less. At these angular scales, fluctuations are mostly 
+3 . dominated by secondary effects due to the local large scale structure (LSS) between us and the last scattering surface 
In addition to generating new anisotropies via scattering and gravitational infall, non- linear effects such as the 
weak gravitational lensing of the CMB Q leave important imprints that can in turn be used to probe cosmology or 
astrophysics related to the evolution and growth of structure (e.g., [||y2|)- 

The lensing effect on the CMB leads to modifications in the direction of photon propagation while leaving the 
surface brightness of the CMB unaffected, implying that its contributions are effectively at the second order level 
in temperature fluctuations. This second order contribution is probed by the three-point correlation function of the 
CMB Q. Here, contributions arise mainly due to the fact that potentials that deflect CMB photons also correlate 
with other secondary effects from the foreground large scale structure, such as the integrated Sachs- Wolfe (ISW; [ jTo| ) 
effect or the Sunyaev-Zel'dovich (SZ; jll]]) effect. At the four-point level, the weak lensing effect itself contributes 
to the correlation function due to its non-linear mode-coupling nature. The trispectrum, the four-point correlation 
function in Fourier space, due to lensing was discussed in |7j and the same trispectrum, under an all-sky formulation, 
is presented in j^] . In ]l3| , we presented the trispectrum due to lensing alone and the correlation between lensing and 
secondary effects and focused on the resulting non-Gaussian covariance of the power spectrum measurements of the 
CMB temperature fluctuations. 

Here, we study an important application of the CMB trispectrum involving a model independent recovery of 
statistics related to the foreground angular deflections of CMB photons (or alternatively the projected potentials that 
generate these deflections). In terms of the trispectrum formulated under a flat-sky approximation, we show explicitly 
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1 http://map.nasa.gsfc.gov 

2 http://astro.estec.esa.nl/Planck/; also, ESA D/SCI(6)3. 
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how one can use a higher order statistic such as the power spectrum of squared temperatures to extract the lensing 
information. The approach described in the present paper was first introduced by Ref. B while other approaches were 
considered in Refs. (UJ^Il based on CMB temperature and polarization observations. We discuss the full trispectrum 
formed by the weak lensing effect on the CMB and extend previous work on this subject by calculating in detail the 
noise associated with lensing extraction. We also discuss noise contributions resulting from the CMB trispectrum 
formed via the correlation between lensing and secondary effects due to the large scale structure. Since the lensing 
effect modifies the damping tail of CMB anisotropics and the information from the damping tail is used to reconstruct 
lensing statistics, we suggest that any secondary effect that dominates small angular scale fluctuations can have a 
significant impact on the lensing extraction. An important contribution at these scales is the SZ effect; fortunately 
this contribution can be removed from thermal CMB maps by noting its distinct spectral signature in observations 
involving many frequencies about the SZ null at ~ 217 GHz |l5| ]. 

As discussed in ]l6| ], the lensing extraction is essential for an important application utilizing CMB polarization 
observations. It is now well known that the curl (or B-modes) of the polarization allows a direct detection of the 
gravitational wave contribution p7j . The amplitude of the angular power spectrum of this polarization signal is 
proportional to the square of the energy scale of inflation, and thus any detection of the B-mode polarization allows 
a probe of the inflationary energy scale. There is one significant source of confusion, however, in the B-mode map 
which can potentially limit the detection of the gravitational wave signal. As discussed in Ref. JlJ] , the weak lensing 
effect on CMB polarization results in a fractional conversion of the dominant E-mode polarization due to density 
perturbations to the B-mode. The contribution from lensing to the B-mode map is such that its power spectrum is 
at the level of the polarization signal expected from gravitational waves if the energy scale of inflation is at scales of 
~ 10 16 GeV. Since we do not know the energy scale of inflation a priori, if it is sufficiently small the lensing effect will 
easily dominate the contribution from primordial gravitational waves. 

In order to identify the gravitational wave signature with no confusion, one must separate the contribution due to 



lensing reliably. As discussed in 
allow a useful separation of the 



L6[ , observations of weak lensing in large scale structure via galaxy ellipticities do not 
ensing contribution. A substantial contribution to the effect comes from potentials 
at redshift greater than 1 to 2, the redshift range currently accessible with galaxy lensing surveys. Even if we have 
a reliable lensing survey out to a redshift of 5, we will only be able to reduce the lensing signal by a factor of ~ 4, 
which is hardly of any use given the order of magnitude change in the gravitational wave contribution that is possible 
if the inflationary energy scale is much less than that related to grand unified theories | JT6| ] . Thus, a better approach 
is to use a source at the last scattering itself, ie. CMB anisotropies, to reconstruct the lensing potential. With 
such a reconstruction, one can clean the polarization maps and obtain unbiased and model independent estimates 
of the polarization signal at the last scattering surface. Here, we discuss in detail the approach related to the 
CMB trispectrum and associated non-Gaussian noises, ignored previously in lensing reconstruction [ pOpl] ] and in 
the application of lensing reconstruction to separate the gravitational wave signature fl6}| . Under certain conditions, 
we find that these additional noise contributions are an order of magnitude below the dominant Gaussian noise 
contribution. If the noise-bias associated with the Gaussian part of the temperature squared power spectrum is 
removed, then the limiting noise level for lensing extraction is determined by the non-Gaussian contribution presented 
here. Even if the noise-bias is not removed, we note that the proposed statistic allows a high signal-to-noise extraction 
of the lensing signal and a separation of the confusing lensed E-mode contribution in the B-mode map from those due 
to inflationary gravitational waves. 

In § ||, we introduce the basic ingredients for the present calculation. The CMB anisotropy trispectra due to weak 



lensing and correlations between weak lensing and secondary effects are derived in § III. In § IV, we introduce a 
statistic involving the power spectrum of squared temperatures, instead of the usual temperature itself, as a probe of 
the angular power spectrum of lensing deflections and exte nd t he discussion to involve polarization in § |yj In § VI 



we discuss our results and conclude with a summary in § VII. The Appendix contains useful formulae relatec 



to 



additional estimators of lensing based on polarization and a combination of temperature and polarization. 

II. GRAVITATIONAL LENSING 

Large-scale structure between us and the last scattering surface deflects CMB photons propagating towards us. Since 
the lensing effect on CMB is essentially a redistribution of photons from large scales to small scales, the resulting 
effect appears only in the second order |l2| . In weak gravitational lensing, the deflection angle on the sky is given by 
the angular gradient of the lensing potential, a(n) = V</>(n), which itself is a projection of the gravitational potential, 



$ (see e.g. |22|), 
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0(m) = -2 I ° dr ^ A }\°, rfar) . (1) 
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The quantities here are the conformal distance or lookback time, from the observer, given by 

and the analogous angular diameter distance 

d A = H^ K 1/2 smh(H n]( 2 r) , (3) 

with the expansion rate for adiabatic CDM cosmological models with a cosmological constant given by 

H 2 = Hi [O m (l + zf + n x (i + zf + n A ] • (4) 

Here, Hq can be written as the inverse Hubble distance today cHq 1 = 2997.9ft.~ 1 Mpc. We follow the conventions 
that in units of the critical density ZH 2 /87rG, the contribution of each component is denoted Vli, i — c for the 
CDM, b for the baryons, A for the cosmological constant. We also define the auxiliary quantities J7 m — Q c + f2& and 
Qk — 1 — 53j ^i) w kich represent the matter density and the contribution of spatial curvature to the expansion rate 
respectively. Note that as f2# —> 0, d A ~> t and we define r(z = oo) = ro- Though we present a general derivation 
of the trispectrum contribution to the covariance, we show results for the currently favorable ACDM cosmology with 
n b = 0.05, fl m = 0.35, n A = 0.65 and h = 0.65. 

The lensing potential in equation (|l|) can be related to the well known convergence generally encountered in 
conventional lensing studies involving galaxy shear |2^| 

= iv 2 0(m) 

r ° , d A {r)d A (r - r) „ _ . . 
dA(ro) 

(5) 

where note that the 2D Laplacian operating on $ is a spatial and not an angular Laplacian. Expanding the lensing 
potential to Fourier moments, 

we can write the usually familiar quantities of convergence and shear components of weak lensing as [fl2| 

7i (A) ± z 72 (n) =-\J ^l^^e^'-^e^ . (7) 

While the two terms k and <f> contain differences with respect to radial and wavenumber weights, we note that these 
differences are artificial and that they cancel under the Limber approximation [p3| . 
The power spectrum of the lensing potential is defined as 

(0(1)0(1')) = (27r) 2 Ml + l')C^, (8) 

where 5p is the Dirac delta function. Expanding the gravitational potential to density perturbations using the Poisson 
equation 24 

and the expansion of a plane wave, we can write the power spectrum of lensing potentials as 

Cf* = - fk 2 dkP(k) [ drxW^fcrJMkn) [ dr 2 W lcn (k,r 2 )ji(kr 2 ) . (10) 
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Here, 

u/Wi. ^ 10 f Ho V G ( r ) d A{r -r) 

and we have introduced the power spectrum of density fluctuations 

(<5(k)5(k')) = (27r) 3 J D (k + k')P(k) , (12) 

where 



2** =^UJ T(fc) ' (13) 

in linear perturbation theory. We use the fitting formulae of [2f| in evaluating the transfer function T(k) for CDM 
models. Here, 8r is the amplitude of present-day density fluctuations at the Hubble scale; with n = 1, we adopt the 
COBE normalization for Sh of 4.2 x 10 -5 pq ], consistent with galaxy cluster abundance (erg = 0.86 ^7|), Note that 
in linear theory, the power spectrum can be scaled in time, P(k, r) = G 2 (r)P(fc, 0), using the the growth function E§] 



G(r) oc 



II (r) ' - ■ " x ' 



Ha J z ( r ) 



d ^ 1 + z '^W)) • ( 14 ) 



In the non-linear regime, one can use prescriptions such as the fitting function by [p9| to calculate the fully non-linear 
density field power spectrum. 

Note that an expression of the type in equation ( [To|) can be evaluated efficiently with the Limber approximation 
p3[. Here, we employ a version based on the completeness relation of spherical Bessel functions [B0[ 



dk^FikMkr^ikr') » Jd- 2 5 D (r-r')F(fc)|, , , (15) 

where the assumption is that F(k) is a slowly- varying function. A similar approach is also presented in Using 
this, we obtain a useful approximation for the power spectrum of lensing potentials as 
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Using equation (^), we note that the power spectrum of convergence is related to that of the potentials through 

C; K = l/4/ 4 Cf . 



III. LENSING CONTRIBUTION TO CMB 

In order to derive weak lensing contributions to the CMB trispectrum, we follow Hu in Ref. |l| and Zaldarriaga in 
Ref. . We formulate the contribution under a flat sky approximation. In general, the flat-sky approach simplifies the 
derivation and computation through the replacement of summations over Wigner symbols through integrals involving 
mode coupling angles. 

As discussed in prior papers weak lensing remaps temperature projected on the sky through the angular 

deflections resulting along the photon path by 

9(n) = 6(n + V0) 

= 6(n) + V^(n)V l 6(n) + iv i 0(n)V J 0(n)V l V J 6(n) + . . . (17) 

As expected for lensing, note that the remapping conserves the surface brightness distribution of CMB. Here, 8(n) is 
the unlensed primary component of CMB, the contribution at the last scattering surface, and 0(n) is the contribution 
affected by the gravitational lensing deflections during the transit. It should be understood that in the presence of 
low redshift contributions to CMB resulting through large scale structure, the temperature fluctuations also include 
a secondary contribution which we denote by s (n) and, in all cases, a noise component denoted by n (n) due to 
instrumental and detector limitations. We denote the total contribution involving these through 9*(n), such that, 
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4 (n) = 8(n) + 8 s (n) + 8 n (n). Since weak lensing deflection angles also trace the large scale structure at low 
redshifts, secondary effects which are first order in density or potential fluctuations also correlate with the lensing 
deflection angle </>. These secondary effects include the ISW and SZ and are discussed in H|l3|]. 
Taking the Fourier transform, as appropriate for a flat-sky, we write 

6(li) = J dhe(n)e- th A 

= &(h)- J |^e(li')£(li,Ii'), (18) 

where 

h') = ^(l! - 1/) (h - h') -h'+\J I^Uai") (is) 

xp(h" + h' - h) {h" ■ h')(h" + h' - h) ■ h' . 

We define the power spectrum and the trispectrum of the 8* fluctuation field, in the flat sky approximation, 
following the usual way 

(Wihj&ih)) =(2^) 2 <5 D (1 1+ 1 2 )C?, 
<8 l (li) • • . 8 i (l 4 )) c = (2 7 r) 2 5 D (l 1 + h + h + U)T*(li, h, h, h) , 

(20) 

where the subscript c denotes the connected part of the above four-point function and the superscript i denotes all 
possibilities including secondary anisotropics and noise. We make the assumption that primordial fluctuations at the 
last scattering surface is Gaussian, such that statistics are fully described by a power spectrum and the non-Gaussian 
four-point correlator, (8(li) . . . 8(i4)) c is zero. A non-Gaussian contribution is generated by lensing effect and other 
secondary contributions. 



(21) 



A. Power spectrum and Trispectrum 

The lensed power spectrum, according to the present formulation, is discussed in |l2| and we can write 

= [i - / (h • i) 2 ] cf + / l^qti^rid - h) ■ hi 2 . 

Thus, the total power spectrum of CMB anisotropies is 

CI = Cf + Cf + Cf , (22) 

where we denote the power spectrum of the noise component by 

Cf - Jsky^V 2 ^ , (23) 

when / s ky is the fraction of sky surveyed, w^ 1 = 47r(s/TcMB) 2 /(£pix^Vpix) |32j is the variance per unit area on the sky 
when ipi x is the time spent on each of iV p i X pixels with detectors of NET s, and df, is the effective beam-width of the 
instrument. The Cf represent the power spectrum associated with contributing secondary effects. In the absence of 
secondary effects, especially in the cases where secondary contributions can be separated from the primary fluctuations 
using statistical signatures or distinct information, such as the spectral dependence of the SZ, we set Cf = 0. Note 
that in general, however, important secondary contributions such as the ISW cannot be separated and will lead to 
additional noise contributions due to correlations with the lensing potentials. 

The calculation of the trispectrum follows similar to the power spectrum. Here, we explicitly show the calculation 
for one term of the trispectrum and add all other terms through permutations. First we consider the cumulants 
involving four temperature terms in Fourier space: 

e(i 1 )...e(i 4 )) c =((e(ii)- J ^e(i 1 ')£(ii ! i 1 ')) (e(h)-J ^e(h')L(h,h'f) e(i 3 )e(i 4 )) + Perm. 

/ (S^ 11 '^ 11 ' 11 ') / ^ ( 1 2 , )i(l2 ! l 2 ')©(l3)©(l4))+Perm. 

(24) 
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As written, to the lowest order, we find that contributions come from the first order term in L given in equation jis|). 
We further simplify to obtain 



e(ix) . . . e(i 4 )) c = ( J pie^^ - h') [(h - 1/) ■ ii'] J i-ke(i a >(i a - h') (h - i a ') • h'Q(h)e(h)) 

= CfCf U{\i + h)<t>(h + hj) [(li + Is) • Is] [(1 2 + 1 4 ) • 1 4 ] + Perm. , (25) 



where there is an additional term through a permutation involving the interchange of li + 13 with li + 14. Introducing 
the power spectrum of lensing potentials, we further simplify to obtain the CMB trispectrum due to gravitational 
lensing: 



T*(h, 1 3 , 1 3 , 1 4 ) = -C£C£ [C^ +h[ (h + h) • h(h + h) • h + C^ +h{ (h + h) • ls(la + I3) ■ Uj + Perm. , (26) 

where the permutations now contain 5 additional terms with the replacement of (£3, 1 4) pair by other combination of 
pairs. 

The non-Gaussian contribution to the trispectrum, through coupling of lensing deflection angle to secondary effects, 
can be calculated with the replacement of, say, 6(13) and ©(I4) in equation ( B4J ) by ©*(l3) and 8'(1 4 ) containing the 
sources of secondary fluctuations, which will produce an addition contribution involving the 6 s contribution associated 
with 6*. Note that we can no longer consider cumulants such as (6(li')6 s (l3)) as the secondary effects are decoupled 
from recombination where primary fluctuations are imprinted. However, contributions come from the correlation 
between © s and the lensing deflection <b. Here, contributions of equal importance come from both the first and second 
order terms in L written in equation (fl9|). First, we note 

(e t (i 1 )...e\u)) c =((eQi)-J |^e(ii')£(ii,ii')) (&(h)-J ^e(i 2 ')L(i 2 ,i 2 ')) e s (i 3 )e s (i 4 )) + Perm 
- -c* (L(\ 2 , -i 1 )e 5 (i 3 )e s (i 4 )) - eg -i 2 )e s (i 3 )e s (i 4 )) 



^j,ef (Z(l 2 , -h')L(h, Ii')e'(l3)e-0U)) + Perm. 



(27) 



Contributions to the trispectrum from the first two terms come through the second order term in L, with the two <j) 
terms coupling to s . In the last term, contributions come from the first order term of L similar to the pure lensing 
contribution to the trispectrum. 

After some straightforward simplifications, we write the connected part of the trispectrum involving lcnsing- 
secondary coupling as 

T et (l 1 ,l 3> l 3> l 4 ) = 

Cf/ Cf {C^(l 3 • li)(l 4 • lx) + Cg(l 3 • 1 2 )(1 4 • la) + [1 3 • (li + la)] [1 4 • (la + h)} Cft+y + [h ■ (h + U)] pa • (la + la)] Cjf 1+U , } 
+ Perm. (28) 

Note that the first two terms come from the first and second term in equation ( p7[ ) , while the last two terms in above 
are from the third term. As before, through permutations, there are five additional terms involving the pairings of 



IV. POWER SPECTRUM OF SQUARED TEMPERATURES 

As written in equation ( p6| ) , the trispectrum formed by lensing alone depends on the power spectrum of deflection 
angle. Thus, with the right configuration for the CMB trispectrum, one can hope to extract some information related 
to lensing, mainly Cf^, given some knowledge on CMB anisotropies, such as Cf. Note that a contribution to the 
CMB trispectrum due to lensing comes from Ct when I = |Ij + lj| when 1^ and lj are the two sides of the quadrilateral 
that form the trispectrum. As illustrated in figure (Q) , this is equivalent to the diagonal of the quadrilateral and can 
be probed easily with a statistic that averages over the configurations for a given length of the diagonal. For this 
purpose, an useful option is to consider the configurations of the trispectrum that contribute to the power spectrum 



G 




FIG. 1. Left: The lensing trispectrum configuration, following equation (|26|), where the diagonal of the configuration contains 
information related to lensing while the sides of the quadrilateral contain information related to primary anisotropy power 
spectra at the last scattering surface. Right: This lensing information can be extracted by a statistic that effectively probes 
the diagonal, such as the power spectrum of squared temperatures discussed here. This statistic sums over all quadrilateral 
configurations, as indicated by thin lines, for a given length of the diagonal. As we find later, the lensing information associated 
with the other — vertical in this diagram — diagonal, proportional to Ct^ +l ,, acts as a source of non-Gaussian noise in 
this extraction. The other contribution to the non-Gaussian noise, through a term proportional to Cii_i 1+ i 2 i, comes from 
interchanging li with 1 — li , which alters the vertical diagonal while keeping the horizontal diagonal fixed. 



of squared temperatures. To see why such a statistic probe the diagonal of the trispectrum, consider the dehnition of 
this squared power spectrum 

(e 2 (l)6 2 (l')> = (2ir) 2 6 D (l + l')C?, (29) 
where the Fourier transform of the squared temperature can be written as a convolution of the temperature transforms 

e 2 (i) = J J^e(ii)e(i-ii). (so) 

Here, it should be understood that 9 2 (1) refers to the Fourier transform of the square of the temperature rather 
than square of the Fourier transform of temperature, [0(1)] 2 . Note that the Fourier transform of the temperature 
fluctuation, in the flat sky, is 

9(1) = J d 2 6e- tl e f{9) . (31) 
To compute the squared power spectrum of CMB temperature, we take 

<e 2 (i)e 2 (i')) = (27r) 2 fe(i + i')c 2 

= /^/ ^<0(ii)©a-iOe(i 2 )e(i'-i 2 )). 

(32) 

Since the lensing power spectrum contributes to the above squared power spectrum via the diagonals of the quadri- 
lateral, the power spectrum of squared temperatures contain terms which are simply proportional to Cf where I is 
the multipole at which the squared power spectrum is measured. As discussed in ||], this then allows an extraction 
of the angular power spectrum of lensing deflections. 

There is also a dominant noise contribution to the squared power spectrum. This noise contribution comes from the 
lowest order correlations and is present even in the presence of Gaussian temperature fluctuations; thus, we denote 
it as the Gaussian noise contribution. In the presence of the non-Gaussian contribution resulting from the lensing 
effect, we write these two contributions as 
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FIG. 2. Optimal filter function W(h + h,h) for the extraction of Cf" from the CMB 2 -CMB 2 power spectrum as a function 
of h and £2- The angle between li and I2 is constrained by the fixed value of I = |li + I2I in each plot. The left plot (I = 100) 
and right plot (I = 2000) correspond to large and small angular scales respectively. The optimal filter removes excess noise at 
multipoles less than ~ 1500 and is weighted to exploit the sensitivity to lensing of the damping tail. 



with 



C? G = 



and the non-Gaussianity captured by the four-point correlator as 

d 2 h r d 2 i 



(2tt) 2 J (2?r) 



2 r(l!, 1 - l Xl l a , -1 - 1 2 ) 



(33) 



(34) 



(35) 



Following our previous discussion on the lensing trispectrum, we can now write the full four-point contribution to 
the squared power spectrum — by replacing li, I2 , 13 and I4 in equation (26) with the required configuration from 
above — such that 



q2NG _ q<, 



+ S W?S (^{^itf+i.i^i • (!i + 1 2)Q C ? + ! 2 ■ Oi + i 2 )Qf ][(i + 1 2 ) ■ (ii + i2)qf +l2l - (h + 1 2 ) • (1 - lOqf^,,,] 
+ c t- i 1+ i 2 | [0 - h + ■ hc% + (1 - h + 1 2 ) • (i - - ii + 1 2 ) • hc% + (i-h + h) ■ (i + h)c® +h] }] . 

(36) 

As written, and discussed above, one of the terms in the non-Gaussian contribution to the squared power spectrum 
is directly proportional to Cf . To extract this angular power spectrum of lensing deflections, consider the approach 
suggested by [|| involving the quadratic statistics formed by the divergence of the temperature-weighted gradients. 
To understand the mechanism involved, consider a filtered version of the squared power spectrum such that 



e 2 (i) 



d 2 h 



TnUi)e(li)6(l - lr) , 



(37) 



where W is a filter that acts on the CMB temperature. This filtering can be either in real space or Fourier, or 
multipolar, space. For simplicity, we have introduced filters that are applicable in Fourier space. To extract the 
optimal filter in the presence of the Gaussian noise alone, we ignore the additional noise contribution from the second 
and third terms of the trispectrum involving integrals over Ct , 1 1 and 



1— li+ial ■ 
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FIG. 3. The extraction of the lensing power spectrum from temperature data alone. We show the dominant Gaussian noise 
(top lines) and the noise associated with the extra terms in the trispectrum due to lensing alone (bottom lines). The curves 
are for three values of the resolution in the CMB temperature map with an effective sensitivity of 1 piK y / sec. 



We can, thus, write the filtered version of the squared power spectrum as 
d 2 h 



rCt q,_ h , W(l, h) [W(-l, - l x ) + W(-l, h - 1)] 



Cf 



J ^ W (I,I 1 )[l-liC? + l.a-li)C|?_ ll |]] [| -^W(-l,-h)[l-hC® + l-(l-h)C®_ h] 



(38) 



We assume that our filter has the symmetry property W(l, li) = W{— 1, — li) = W{— 1, li — 1) and then check to see 
if this assumption is justified. To be consistent with prior results, we also make one more assumption here. We 
take the Gaussian noise to be the dominant noise contribution and ignore the noise from the two extra terms in the 
non-Gaussian trispectrum that involves integrals over Ct^ + y, and CjJ_ li+ | Thus, the filter we will derive here is 
only optimal in the presence of the dominant Gaussian noise and is sub-optimal when all noise terms are considered. 
We will later add in these two terms as an extra source of noise, with the filter derived by ignoring them applied 
appropriately. 

With these assumptions, we can write the signal-to-noise ratio for the detection of the 4> power spectrum per each 
multipole moment as, 



(21 + 1) 



\ r 4 



/ ©rtm, ii)[i • iicg + 1 • (i - ioq?_ U |] + / m*w& h) 2 2Cfc^ hl 



(39) 



The optimal filter in the presence of Gaussian noise is the one that will maximize this ratio. 

To extract the required filter, we take derivatives with respect to W and find that the maximum is satisfied when 



W(l,h) 



[i-i a c® + i- ,] 



(40) 



Performing the explicit substitutions 1 — > —1, li — > — li and 1 — > —I, li — »• li — 1, we see that our optimal filter does 
possess the desired symmetry. We introduced this filter in Ref. when discussing the CMB 2 -SZ power spectrum as a 
way to extract the lensing-SZ correlation directly using the power spectrum of squared temperature and temperature. 
It is not a surprise that the same filter appears now when extracting the lensing-lcnsing correlations using the power 
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spectrum of squared temperature and squared temperature. A physical mechanism to achieve this filter in real space 
is presented in Ref. |2p[] by taking the gradient of the temperature, weighting it by a noise filtered CMB map, then 
taking the divergence of this product. 

With the filter applied, we get a direct estimator of C, . Introducing, 



Nr 1 



d 2 h [MxC® +1.(1-10^, 



1 ~j (2-) 2 nqui 

(41) 

The filtered squared power spectrum is such that 

N? Cf =Ni + Cf* , (42) 

where the Gaussian noise associated with the extracted Cj is A/. 

We can now add in the noise contribution arising from the two extra terms in the non-Gaussian trispectrum which 
we dropped when deriving the optimal filter in the presence of the Gaussian noise only. This extra noise contribution, 
neglected in prior studies, is 

{cff +h[ [h ■ (h + h)c® + 1 2 • (ii + i 2 )cg][(i + 1 2 ) • (h + i 2 )cjf +Ia | - (h + h) • (i - ioqt,,,] 

+ c |t-i 1+ i 2 | K 1 ~h + h) ■ hCf 2 + h) ■ (1 - IiJCf-mlhO - Ii + la) ■ liC? + (1 - li + la) • (1 + h)Cf l+h] ]) . 

(43) 

With the non-Gaussian noise, we can write 

NfCf = Ni + A, NG + C+* . (44) 

There is a strong possibility that the noise-bias Ni associated with the Gaussian contribution to the squared power 
spectrum can be removed during the lensing reconstruction. This would involve the construction of a new estimator 
of lensing, or filtering scheme, that removes iVj and leaves iV ; NG as the final noise contribution associated with the 
reconstructed Cf^. The removal of the dominant Gaussian noise bias is extremely useful since this allows one to lower 

the noise contribution associated with Cf by at least an order of magnitude. We will return later to the impact of 
such a noise-bias removal on the identification of an inflationary gravitational wave signature in the CMB in the form 
of a curl (or B-mode) contribution to the polarization. 

In addition to the additional non-Gaussian noise contribution from the lensing trispectrum itself, there is an 
additional noise contribution resulting from the relevant configuration of the trispectrum due to lensing-secondary 
correlations. We do not reproduce this contribution here, but it can be easily calculated by replacing (li, I2, 13, 14) in 
equation (p8[), with (li, 1— li, I2, — 1— 12). Note that this additional noise results from the correlation between potentials 
that deflect CMB photons and the secondary contributions associated with the ISW and SZ thermal effects. In the 
presence of such contributions, we also note that Cf in equation ([ll]), contains all contributions to noise, including 
the secondary contribution. It is well known that certain secondary contributions can be removed based on distinct 
signatures, such as the spectral dependence of the SZ effect Jl^]. In such cases, we expect only the lensing-ISW 
correlation to produce a noise contribution. This effect is small and can be easily ignored. In the presence of the SZ 
contribution, however, there is a significant noise contribution that can potentially limit the extraction of the lensing 
power spectrum. 

Note that the filter is symmetric in li and 1 — li. We show surfaces of this filter in figure (||) as a function of l\ 
and I2 = |1 — li| for two values of I corresponding to large and small angular scales respectively. As shown, the filter 
behaves such that information related to lensing is extracted at arcminute scales via the damping tail of the CMB 
anisotropics. Additionally, since lensing smooths the acoustic peak structure, further information is also extracted in 
a subtle way from fluctuations associated with large angular scales. Thus, a simple filter which effectively cuts off the 
acoustic peak structure but extracts information from small angular scales, such a a step function in multipolar space, 
does poorly when compared to the filtering scheme proposed above p3[ . As mentioned before deriving the filter, it 
should be noted that the proposed filter is only optimal to the extent that Gaussian noise contribution is dominant. In 
the presence of noise contributions from additional four-point correlations or, more importantly, additional terms in 
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FIG. 4. The extraction of lensing power spectrum from temperature data in the presence of foregrounds (mainly the thermal 
SZ effect). The increase in dominant Gaussian noise (top lines) is due to the extra SZ power spectrum, which peaks at small 
angular scales. The dot-dashed line is the resulting non-Gaussian noise contribution due to the trispectrum formed by the 
SZ-lensing cross-correlation. 



the lensing trispectrum itself, the filter is sub-optimal. We failed to find a simple analytical description for an optimal 
filter that also suppresses the non-Gaussian noise contribution from the additional terms in the lensing trispectrum. 
As we find later, since these terms only contribute at an order of magnitude below the dominant Gaussian noise, we 
do not expect this to be a significant problem in the extraction of the lensing information. 



V. THE EB TRISPECTRUM 



In addition to temperature, lensing also affects the CMB polarization. We discuss here the best estimator from 
polarization observables for extracting the lensing contribution following the discussion in Ref. pH and provide 
relevant equations for other estimators in the Appendix. The statistic considered in this section involves the squared 
map constructed by convolving, in Fourier space, the E- and B-mode maps of the polarization. We denote this by 
(EB)(l). To understand this combination, we consider the Fourier space analogue and, as before, write, 

EB® = \ J h) + B(h)E(l h)] . (45) 

Now considering the power spectrum formed by (EB(\)EB(V)) = (2tt) 2 (5(1 + V)C\ ' , we write 

{EB{\)EB{\')) = lj^J f^L {mi )B{\ h) + B(h)E(l - h)\ [E{h)B{Y 1 2 ) + B{h)E{\' 1 2 )]) . (46) 

Here, again, we find two sources of contributions, involving a noise contribution due to Gaussian fluctuations in the 
E- and B-maps, 

r ,{EB) 2 G _ 1 f d 2 li \ r>r ,EEtriBBt , onBBtrtEEt 1 / A n\ 

; ~~ 4 J (2^r)2 L h I 1 " 11 ' ll MilJ ' ^ ' 

and a non-Gaussian contribution due to the lensing effect on polarization. 

To calculate the lensing contribution to the polarization, we follow the notation in Ref. Introducing, 



±X(n) = ±X(h + V<£) (48) 

2 



±X(n) + V^(n)V l ±A(n) + -V^(n)V^(n)V l V J ±A(n) , 
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FIG. 5. The extraction of lensing power spectrum from polarization data. We show the associated noise due to the dominant 
Gaussian noise (top lines) and noise due to extra terms in the trispectrum (bottom lines). The curves are for a resolution of 5 
arcmin in the CMB temperature map with an effective sensitivity of 1 fj,K y^sec. As shown, the EB estimator probes the lensing 
potential power spectrum to smaller angular scales than the quadratic estimator on temperature data alone or estimators based 
on other combinations of polarization and the temperature. 

where ±X = Q±iU represent combinations of the Stokes parameters. The lensing effect is again to move the photon 
propagation directions on the sky. In Fourier space, we can consider the E- and B-mode decomposition introduced in 
Rcf. [O such that ±X(1) = 25(1) ± iB(l) and 



1(1) = ± X(1) - J ^ ±X (h)e ±2i ^-^L(l, h) 



(49) 



With this, we can write 
25(1) = 25(1) - 
5(1) = B(l) - 



d 2 h 



(2tt) 2 



[25(li) cob2(( P i 1 - tpi) - B(h) sin 2{tp h - <&)] 0(1 - l x )(l - h) ■ h 
2 [25(h) sin 2^ -<pi) + B(h) cos 2^ -(pi)] 0(1 -li)(l- U) - U , 



(50) 



Since primordial 23(1), due to the gravitational wave contribution, is small, we can make the useful approximation 
that 



25(1) = 25(1) 



d 2 h 

(2^) 2 



[25(h) 0082^ - <&)] 0(1 - lr)(l - h) • h 



2?(1) = - 



d 2 h 

(2^) 2 



[25(h) 8^2(^1, - (pi)] 0(1 - li)(l - li) ■ li , 



(51) 



Under this approximation, the lensed polarization power spectra can be expressed in terms of Cf and the unlensed 
quantities following Ref. [O : 



l 



i r d% 



Ct + 2 J (^2C | r ii| [(l-l 1 )-l 1 ]^Cr[l + cos(4^ 1 )], 



Cf B = \ J ^C,f h] [(l - h) ■ k] a Cg*[l - cos(4^)] 



(52) 



The total polarization power spectra is the sum of these lensed spectra and a noise contribution 
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BBt 



c, 



BB 



Cf, 



(53) 



where the noise component is given by equation ( |23| ) with a possibly different NET s for polarization. 

With the lensing effect on -B(l), we can now calculate the trispectrum contribution to the squared power spectrum 

as 



(BB) NG 



d 2 h 

(2tt)2 



1 r 4> 

d% rl 



( 2;r) 2 { 4"l 1 i+ 1 ^l 



c^c^ih + h) ■ (i + i 2 )(i! + 1 2 ) ■ h - c^c$x { (h + h) ■ h(h + h) ■ (i - h) 



x sin 2(^1^ - <pi+i 2 ) sin2(^ h + ip h ) 



4 °|i-ii+i 2 | 



c^Xf^ M {\ - ii + h) • (i + i 2 )(i - ii + 1 2 ) ■ (i - h) - c?»c?»(i -h + h) ■ h(i -h + h) ■ h 

x sin2((^i 1 - < y 5i + i 2 )sin2(v3i_i 1 + (pi 2 )j , 



(54) 



To extract the lensing contribution, one can again consider a filtering scheme, and following arguments similar to 
the case with temperature, the derived filter that maximizes the signal-to-noise in the presence of Gaussian noise is 



C?v\ ■ h + I ■ (1 - l x ) sin 2^ + <p^ h ) 



°h °ll-lll 



(55) 



With this filter applied, the associated noise due to the dominant Gaussian part is Ni, where, 



ivr 



d 2 h 



CfH ■ h + C*X { 1 • (1 - l x ) sin 2 2(tp h + w- h ) 



(2nf 



r,r<EEtr'BBt I r >r <BBtr<EEt 

zu h °ll-lll + ZLy h °ll-lll 



(56) 



The extra noise due to the additional terms in the trispectrum is 



1 



4 lh+12 



(2tt) 2 J (2tt) 

c^cf^ih + h) • (1 + l 2 )(li + 1 2 ) ■ li - C^Cf^ h] (h + h) ■ h(h + h) ■ (1 - li) 



x sin2(y)i_i 1 - ^i + i 3 )sin2(<^i 1 + ip h 



+ 



4 w |i-ii+i 2 



Cf?X\ C \L+h\( l - *i + ^ ' + W - li + la) • (1 - li) - CT'CT(1 - li + I2) • li(l - li + h) ■ I 



xsin2(<^i 1 - </?i+i 2 ) sin2(<^i_i 1 +<^i 2 )| 



(57) 



VI. DISCUSSION 



We have introduced the trispectrum, or four-point correlation function, of CMB anisotropics as a probe of the 
large scale structure weak lensing effect. The trispectrum in CMB anisotropics is generated by the non-linear mode 
coupling behavior of the weak lensing effect on the CMB. As written in equation (p6|), the trispectrum formed by 
lensing contains a significant contribution from Cf when I = \\\ + 13 1 . li and I3 are the two sides of the quadrilateral 
that formed the trispectrum in the Fourier space as illustrated in figure |l|. Contributions thus probe the lensing 
effect via the diagonal of the trispectrum configuration, and statistics sensitive to the diagonal will provide a practical 
approach for extracting the lensing angular power spectrum. 

A useful statistic for this purpose is the power spectrum formed by squared temperatures introduced earlier. This 
approach expands the suggestion in Ref. [p3fl to use the squared temperature-temperature power spectrum to study 
the CMB bispectrum and a filtered version to extract the power spectrum associated with lensing-SZ correlation. The 
power spectrum formed by squared temperatures probes the trispectrum. This power spectrum, however, contains 
a dominant contribution due to the Gaussian fluctuations associated with temperature which potentially dominates 
the contribution arising from the non-Gaussian weak lensing effect. To optimally extract the lensing contribution 
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FIG. 6. Left: Weak lensing as a CMB experiment. We show errors for the reconstruction of convergence, or projected 
mass density, to the last scattering surface with Planck temperature data. This is compared to the convergence measurements 
expected from observations of galaxy ellipticities which probe to redshifts of ~ 1 to 2. Right: The cumulative signal-to-noise 
ratio for the detection of the lensing effect in Planck data, using the filtered squared temperature statistic discussed here (solid 
line) and the gradient statistic (dotted line) described in Ref. [14] . The dot-dashed line is the cumulative signal-to-noise ratio 
when the dominant Gaussian noise is removed as part of the lensing reconstruction. With such a removal, the noise in the 
lensing reconstruction is determined by the non-Gaussian contributions which are roughly an order of magnitude below the 
Gaussian noise level. 



in the presence of this dominant Gaussian noise, we have introduced a filter in Fourier space. This filter requires 
knowledge of the noise contributions to the map as well as the power spectra of prelensed primary anisotropies at the 
last scattering surface. Though the filter can be easily applied to the Fourier transform of the squared temperature 
map in the Fourier space, a useful direct approach is presented in Ref. p0| . This real space construction involves 
the divergence of the temperature-weighted temperature gradients, with each of these temperature estimators first 
filtered for the excess noise. 

We have extended the previous discussions of the extraction of lensing signal to include a discussion of the full 
trispectrum. The full trispectrum contains terms which were previously ignored and these terms contribute as addi- 
tional sources of noise in the lensing reconstruction and can potentially interfere with the extraction of the lensing 
signal from CMB data. These additional non-Gaussian noise terms correspond to the lensing contribution associated 
with the second diagonal of the trispectrum configuration shown in figure (|l]) and the same diagonal under a permu- 
tation, while keeping the same vector 1 for the diagonal from which lensing information is extracted. In figure (^), 
we summarize our results related to the lensing extraction using temperature data alone. We show the dominant 
Gaussian noise contribution and the associated non-Gaussian noise contribution, due to the two additional terms in 
the trispectrum which involve integrals over Ct^ +l , and Ct_ li+l ,. We show these noise contributions for temperature 

maps constructed with an effective detector sensitivity of 1 fjK y / iec and varying resolution with beam widths of 5 to 
1 arcmin. As shown, non-Gaussian noise contributions from additional terms in the trispectrum are roughly an order 
of magnitude below the dominant Gaussian noise contribution. 

This conclusion is subject to several conditions. As summarized in figure (0), the extraction from temperature 
alone can be significantly impaired by the presence of foregrounds that correlate with the lensing deflection angle. For 
example, the SZ contribution hinders the separation in two ways. The first is an enhancement of the Gaussian noise 
signal by the addition of Cf z to C\ as indicated by the difference between top dashed and the dot-dashed lines. The 
second is a noise contribution due to the non-Gaussian trispectrum formed by the correlation of lensing potentials 
with the SZ effect due to the large scale structure. We show this with the long-dashed line. In the presence of SZ, 
it is now clear that the lensing extraction is strongly limited. The other source of correlation is the ISW effect. Its 
effects, through the contribution to Gaussian noise, are already included in the Cf defined earlier. The power spectra 
calculated using the publicly available CMBFAST code |34| contain the ISW signal as well. The noise due to the 
correlation between ISW and lensing is significantly smaller and is at a level below 10 -10 in figure (Q). Thus, one can 
easily ignore the additional noise contributions resulting from the ISW effect. Furthermore, the SZ noise contributions 
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can also be controlled. One can easily separate the SZ contribution from thermal CMB fluctuations based on the 
frequency spectrum of the SZ effect |L5|]. In the upcoming Planck mission, SZ separation based on multifrequency 
data will prevent this contribution from adding significant noise to CMB maps. We suggest that small angular scales 
experiments seeking to extract the lensing signal based on temperature anisotropies be designed to contain channels 
that are optimized for a separation of the SZ signal from thermal CMB fluctuations. 

In addition to temperature, lensing also modifies polarization data. Thus, one can employ estimators based on 
polarization, and combinations of polarization and temperature, to extract information related to lensing. In figure (|^), 
we illustrate the usefulness of polarization observations for a lensing separation. The dotted lines are the errors 
associated with the extraction based on the EB estimator. This estimator was suggested by ]2l| as an estimator that 
probes the angular power spectrum of lensing on scales smaller than that probed by temperature or other combinations 
of temperature and polarization. EE and 0E estimators also probe the lensing deflection power spectrum, and 
have Gaussian and non-Gaussian noise contribution in between those of the 08 and EB estimators. While we 
presented a detailed derivation of the EB estimator, we write down the trispectrum and related filters for other 
estimators associated with polarization and combinations of polarization and temperature in the Appendix. As 
shown in figure (^|), we again find the non-Gaussian noise contributions due to additional noise terms in various 
trispectra due to lensing alone to be roughly an order of magnitude below the dominant Gaussian noise contribution. 

Similar to the 00 estimator, the OE estimator contains an additional noise term due to the correlation between 
weak lensing and secondary anisotropies such as the SZ effect. We have presented a derivation of this contribution in 
the Appendix. In both these cases, note that we have only considered the secondary contributions to the temperature. 
Unlike contributions to the temperature, secondary contributions to the polarization are significantly smaller and are 
higher order in density fluctuations; for example, in addition to the density field, these secondary effects also depend 
on the velocity field [ [35[ . This makes the secondary polarization-lensing correlations less significant, similar to the 
correlation between lensing potentials and the moving baryons responsible for the Ostriker-Vishniac/kinetic SZ effect 
p6| . Therefore, we ignore noise contributions from secondary effects associated with lensing extraction based on 
polarization estimators alone. 

In general, the CMB data allow high signal-to-noise ratio reconstruction of the lensing effect. We illustrate the 
extent to which Planck's frequency cleaned temperature maps can be used as a lensing experiment in figure (^). Here, 
we show the angular power spectrum of convergence for sources at the last scattering, ie., CMB, and at a redshift 
of 1; the latter is consistent with what is expected from weak lensing surveys via galaxy ellipticity data. In the 
right panel, we show the cumulative signal-to-noise ratio for Planck based on the quadratic statistic discussed here 
and the gradient-based statistic of 14 1. The reconstruction based on the squared power spectrum has factor of ~ 3 
higher cumulative signal-to-noise than the one based on gradients in the temperature map. This can be understood 
by noting that the latter method only extracts lensing from the first order expansion of the temperature in terms of 
the deflection angle, while significant contributions also come from the second order term that is also probed by the 
proposed estimators. 

As shown in figure (^|), due to Planck's beam of ~ 5 to 10 arcmin the reconstruction of convergence becomes noise 
dominated at multipoles of ~ 1000, while one can probe lensing to much smaller angular scales with galaxy information. 
The difference between the two curves illustrate the importance of lensing reconstruction for the polarization studies 
related to gravitational waves. Since CMB photons are affected by potentials out to the last scattering surface, the 
resulting lensing effect is at least an order of magnitude higher than the lensing effect on sources that are visible from 
the large scale structure out to a redshift of a few. Thus, it is unlikely that observations of the large scale structure 
lensing can be used to subtract the lensing contribution to the B-modc polarization sufficiently to permit a confusion 
free detection of the polarization signal due to gravitational waves. We refer the reader to |ll| for a further discussion 
on this application. 

We can extend the discussion in Jl6| ] to consider one interesting possibility. As discussed earlier, one can make a 
significant improvement in the lensing reconstruction when the Gaussian noise bias associated with the temperature 
squared power spectrum is removed. In such a scenario, the noise contribution associated with the lensing extraction 
is determined by the non-Gaussianities which are usually an order of magnitude below the dominant Gaussian noise 
contribution. Such a Gaussian noise bias removal allows an extremely high signal-to-noise extraction of the lensing 
signal even for experiments with poor sensitivities that are not optimized for the lensing reconstruction. As shown 
in figure (|^) , for Planck, the removal of the dominant Gaussian noise leads to detection of the lensing angular power 
spectrum with a cumulative signal-to-noise ratio of ~ 700, which is roughly a factor of 20 above the level expected when 
the dominant Gaussian noise contribution is not removed. Such an extraction not only improves the effectiveness of 
cosmological studies involving lensing of the CMB jj7| , but also allows a high signal-to-noise separation of the lensed 
E-mode contribution to the B-mode polarization map of the CMB. This improves the confusion- free limit on the 
inflationary energy scale attainable with CMB polarization experiments by another order of magnitude |l6[ ] . 

Following , we illustrate the importance of noise-bias removal in figure (0). Here, we compare various contribu- 
tions to the power spectra of CMB polarization, including the lensed E-mode contribution to the B-mode polarization 
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FIG. 7. Left: Various contributions to the CMB polarization. The long-dashed curve shows the dominant polarization signal 
in the gradient, or E-mode, component due to scalar perturbations while the solid line shows the curl, or B-mode, polarization 
signal from the gravitational-wave background assuming an inflationary energy scale of 2.3 x 10 16 GeV, consistent with the 
limit based on COBE. The short-dashed curve shows the contribution to the B-mode power spectrum due to gravitational 
lensing. The dotted curve labeled 'G' is the residual on this contribution due to Gaussian noise when lensing deflections are 
estimated from a temperature map with a resolution of 1 arcmin, a sensitivity of 1 fiK -y/sec and a sky coverage of 80%. The 
dot-dashed curves labeled 'NG' is the residual if the dominant Gaussian noise is removed as part of this extraction leaving only 
the non-Gaussian noise of equation ( pl3| ) . Right: Minimum inflation potential observable at la as a function of survey width 
for a one- year experiment. The solid curve shows results assuming no lensed E-mode contribution to the B-mode map while 
the short-dashed curve shows results including the effects of lensing as a noise contribution. The dotted and dot-dashed curves 
correspond to the lensing subtraction described in the left panel. The order of magnitude improvement in the inflationary 
energy scales that can be probed confusion- free with CMB polarization data suggests the importance of investigating ways to 
remove the noise-bias associated with the lensing extraction based on the squared power spectrum. 



map. We also show the residual lensing contribution to the B-mode map after the lensing effect has been largely 
removed using an estimator for lensing deflections. The estimator employs a CMB temperature map with a resolution 
of 1 arcmin and a sensitivity of 1 /xK y/sec. The curve labeled NG is the residual contribution when the lensing extrac- 
tion is carried out to the noise level determined by the non-Gaussian noise contribution, assuming that the dominant 
Gaussian noise contribution can be removed during the lensing extraction. As shown, such a removal provides an 
order of magnitude improvement over the residual from the dominant Gaussian contribution. 

In the right panel of figure (Q), we show limits on the inflationary energy scale based on lensing extraction with 
the remaining residuals depicted in the left panel of the same figure. We refer the reader to jl6| for details related 
to this calculation. For comparison, we also show the inflationary energy scale limit when the lensing contribution is 
ignored (only instrumental noise) and the limit when the lensing contribution is treated as an additional contribution 
to the noise associated with the B-mode map. As illustrated in these figures, the potential removal of the Gaussian 
noise associated with the lensing reconstruction reduces by an order of magnitude the confusion-free upper bound 
on the gravitational wave signal in the B-mode polarization. Such an improvement motivates the development of 
practical schemes to remove the noise-bias associated with the dominant Gaussian part of the temperature squared 
power spectrum. In future work, we will return to this issue and will investigate ways to achieve the required level of 
noise removal. 
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VII. SUMMARY 



We have discussed the four-point correlation function, or the trispectrum in Fourier space, of CMB temperature 
and polarization anisotropics due to the weak gravitational lensing effect by intervening large scale structure. The 
trispectrum contains information related to the lensing deflection angle, which for a given configuration of the trispec- 
trum is captured primarily by the diagonal of the quadrilateral that defines the trispectrum in Fourier space. We 
have devised a statistic involving the power spectrum of squared temperatures that captures this information from 
the diagonal. This provides an observational approach to extracting the power spectrum of deflection angles induced 
by the weak gravitational lensing effect on the CMB. The method extends previous suggestions on the use of the 
squared temperature power spectrum to extract the lensing-secondary correlation power spectrum by probing the 
bispectrum formed by the lensing effect. The technique presented here is not new and has been discussed previously 
in the literature by p0[ . We extend this work by providing a full derivation of the trispectrum terms both due to 
lensing and lensing-secondary correlations, and by providing a detailed discussion of the filtering required to mini- 
mize the dominant noise component due to Gaussian temperature fluctuations that impairs the extraction of lensing 
information. 

We also extend previous discussions on the trispectrum and associated weak lensing reconstruction from CMB 
data by calculating in detail contributions to the noise in the reconstruction. Previous work discussed the dominant 
Gaussian noise, but we include in addition terms associated with non-Gaussian contributions due to lensing alone and 
the correlation of the lensing effect with other foreground secondary anisotropies in the CMB such as the ISW effect and 
the SZ effect. When the SZ effect is removed from the temperature maps using its spectral dependence, we find these 
additional non-Gaussian noise contributions to be an order of magnitude lower than the dominant Gaussian noise. 
When SZ is not removed, the extraction of lensing information is significantly suppressed; we strongly advise the use of 
multiple frequencies in upcoming small angular scale CMB experiments, if there is an interest in using data from such 
experiments to conduct a weak lensing study. The lensing extraction provides a high signal-to-noise ratio detection of 
the lensing deflection power spectrum. For the Planck surveyor, with frequency-cleaned t emp erature maps, detection 
will have a cumulative signal-to-noise of ~ 35, while certain techniques in the literature [ft4|, involving the gradients 
of the temperature map, only give a cumulative signal-to-noise ratio of ~ 10. Since the extraction with Planck is 
severely limited by its beam of 5 to 10 arcmin, we strongly recommend higher resolution observations at the arcminutc 
scale. If the noise-bias associated with the dominant Gaussian part of the temperature squared power spectrum is 
removed, Planck can extract the lensing power spectrum with temperature data alone with a cumulative signal-to- 
noise ratio around ~ 700. With information from polarization, the signal-to-noise ratio associated with lensing can 
be further increased by another factor of a few. The order of magnitude improvement in the lensing separation 
motivates the development of practical methods that can remove the dominant Gaussian noise contribution to the 
squared temperature power spectrum. In general, our study suggests that future CMB observations can be used as a 
weak lensing experiment akin to current attempts involving lensing reconstruction from galaxy ellipticity data. This 
opportunity will help to achieve an ultimate challenge for CMB experiments, a confusion-free direct detection of the 
gravitational wave signature from inflation. 
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VIII. APPENDIX I. EE AND 6E ESTIMATORS 

Here we provide the appropriate formulae for deriving optimal filters for the EE and OF, estimators of Cf^ . 

A. EE estimator 

The Fourier transform of the squared E-mode map can be expressed as a convolution of Fourier transforms of the 



E-mode map itself, analogous to equation ( p0|) of § [IV 
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(EE)(l) = J -j0-E(h)E(l-h) 



(58) 



Defining the power spectrum for the squared E-mode map as (EE(Y)EE(Y)) = (27r) 2 5(l + Y)C E E , we write 
(EEQ)EE(Y)) = J^2 J -^{E(h)E(l-h)E(h)E(l r -1 2 )) 



(59) 



As before, this power spectrum will consist of a noise component due to Gaussian fluctuations in the E-mode map 
and a non-Gaussian contribution due to lensing given respectively by, 



E Z E Z G 
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x cos 2(^1! - Lp\+\ 2 ) cos 2(<^i_i 1 + tpi a ) I . 



(61) 



In an analysis similar to that of § IV and § 0, we find that the optimal filter in the presence of the Gaussian noise 
term above is 



W(l,h 



c EE \ ■ h + qf4,i • (1 - Ix) cos2( Wl + w-h) 



onEEtr'EEt 
zu h °|i-ii| 



Using this filter, the dominant Gaussian noise Ni is given by 



(62) 



d 2 h 



C EE \ ■ h + C EE J ■ (1 - h) cos 2 2( Wl + tpi- h ) 



(2n) 2 



onEEtriEEt 
ZL, h °|l-li| 



while the excess noise due to the additional terms in the E-mode trispectrum is 

fh_ t d 2 l 

(2^ J (>) 2 



CZ+xAC^h + C^-ih + h) c^ E 2l (i + h)-c EE hl (i-h) -Oi + l 



x cos 2(y>i_i 1 - (pi +h ) cos 2((p h + ip h ) 



c; 



ll-ll+l 



C|f4,(l - h) + C EE h] • (1 - h + h) [C,f4|(l + l a ) - C EE h •(!-!! + 1 2 ) 
x 0082(^1! - <^i + i 2 )cos2(^i_i 1 + ipi 2 )\ . 



(63) 



(64) 



B. OE estimator 



Here we derive the optimal filter and accompanying noise for the OE estimator of Cf . The following equations ( pq ) 
through ( |7l|) for this estimator are entirely analogous to equations (|5^) through (^4|) of the preceding subsection of 
this Appendix. An additional complication here is that, unlike in the EB estimator where E and B are uncorrelated 
due to parity, 9 and E are correlated with a power spectrum defined by Cf E below. This produces an additional term 
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through the lensing effect in addition to that involving correlations between 9 and E separately. First, we expand in 
Fourier space the product of 9 and E maps, 



(65) 



to obtain 



(6^(1)6^(1')) = J / 7^2 / J^><[ e (W ^) + Ml Pfo W la) + £(1 2 )6(1' - 1 2 )]) , (66) 



and 



(2nf 



/0 h °|l-li| + ZL/ ii °|l-li|+ 4 °h ° 1 li 



(67) 



The non-Gaussian part is 



C, 



(6E) 2 NG 



+ l/(Sp/ K 1 * -(h+h) [qff l2| (l + l 2 )-C,f_ s li| (l-l 1 )].(l 1 +l 2 )cos2(^_ ll -^ +l2 ) 

+ [C E % + C EE \ 2 ] ■ (h + h) [q? +l2 ,(l + h) C®_ hl (\ h)] ■ (h + h) cos2( Vll + Vl2 ) 

+ [Cf E cos 2^ + ^ 2 )l! + C% E h] ■ (h + h) [qf4| cos2(^,_ ll - mi2 )(l + 1 2 ) - 0^,(1 - li)] • (li + 1 2 ) 

+ [C£ £ li + Cg E cos2(^ ll + ^ l2 )l 2 ] • (h + 1 2 ) [e®f l2| (l + 1 2 ) - qf^i cos2(^ 1 _ ll - ^ 1+ , 2 )(1 - h)] ■ (h + 1 2 )} 

+ C \A+h\{ [ C ®-h\( l - li) + cos2(^ 1 _ ll + ^, 2 )1 2 ] ■(!-!!+ 1 2 ) [q?f la| (l + 1 2 ) - Cf E cos2(^ ll - m , 2 )l 1 ] ■ (1 - h + 1 2 ) 

+ [Cjf^i cos2(^ 1 _ ll + Vla )(l - lx) + C£ B 1 2 ] ■ (1 - h + h) [e^f l2| cos2(^ ll - mi2 )(l + 1 2 ) - cf E h] ■ (1 - h + h) 

+ [c EE hl (l + C EE h] ■ (1 - li + h) [C|?+i 2 |(l + 1 2 ) - Cffh] ■ (1 - li + 1 2 ) cos 2^ + Vla ) 

+ [qt^i (1 -h) + Cf\ 2 ] ■(!-!!+ 1 2 ) [e,f + s l2| (1 + 1 2 ) - C EE h] ■ (1 - U + 1 2 ) cos2( m , 2 - ^ij} } . 



(68) 



The filter required to suppress the dominant Gaussian noise and extract lensing information is 

'cf E \ ■ h + C® E h \ l ■ - li)l [1 + cos 2^ + 



W(l,h) 



and the associated Gaussian noise is 



r<etr<EEt I riEEtr<Qt i or-efit^efit 



(69) 



d 2 h 



Cf E \ • li + e.ff,,! • (1 - Ix) [1 + cos 2^ + ^_ h ) 



(2tt)2 



ZO Zi °|l-li| + °|l-li| +4 °ii °|l-li| 



(70) 



while the non-Gaussian noise is 



x { C |t?+i 2 |{ K 1 ^ + C ?M ■ (!i + h) [C,f_fi 2 |(l + 1 2 ) - e^^l - h)] • (h + 1 2 ) cos2(^ 1 _ ll - <p l+h ) 
+ [C EE h + C EE h]-(h+h)[c® +hl (l + h)-C®^ 

+ [e^cos2(^ ll + ^ l2 )i 1 + cf E \ 2 ] • (h + h) \c® E hl cos2(^ 1 _ ll - mi2 )(i + 1 2 ) - qfi^i - i a )l • (h + 1 2 ) 
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+ [C£ B li + Cf E cos2( v3ll + ^ 2 )1 2 ] • (li + 1 2 ) [qtf l2| (l + 1 2 ) - Cjf^i cos2(^ 1 _ ll - mi2 )(l - h)] ■ (h + 1 2 )} 
+ C \i-i 1+h \{ [q?4|(l-li) + C® E cos2( W _ ll+ ^ 2 )l 2 ] .(1-lx + la) 



qT+i 2 |(l + la) - cos 2 (^ii - W+OUJ ■ (1- li + la) 

cos 2^ + ^ l2 )(l - h) + Qf E l 2 l • (1 - li + l a ) [cff+L cos 2^ - ^ 1+l2 )(l + 1 2 ) - C° E h] ■ (1- h + 1 2 



Cji_i x |(l - li) + CE*l a j • (1 - li + l a ) [q? +Ia |(l + la) - C^liJ • (1 - li + la) cos 2^ + ^ 2 ) 
Cft^l - h) + Cf 2 h] •(!-!! + 1 2 ) [qf+ £ l2 |(l + 1 2 ) - C E %] ■ (1 - la + l 2 )cos2( m , 2 - Wl )}} . 



The trispectrum formed from two temperature and two E-mode terms will contain an additional non-Gaussian 
contribution through the coupling of the lensing deflection angle to secondary effects. A similar contribution occurs in 
the temperature trispectrum and was discussed at the end of § [II. As in equation (p7|) , cumulants such as (0(1)0 S (1')) 
will vanish as the secondary effects are decoupled from the primary fluctuations generated at the surface of last scatter. 
However, contributions come from the correlation between s and the lensing deflection <f>. As before, contributions 
of equal importance come from both the first and second order terms in L written in equation (^). We consider this 
additional non-Gaussian contribution for the connected portion of a four-point function in Fourier space 



(71) 



^(i 1 )^(i 2 )e t (i 3 )e t (i4)) =( £(10- 



(2?r) ; 

721/ 



■S(l 1 ')cos2( Wi - tp h )L(h,h r 



E(h)~ J ^ J E(l2')cos2( W3 - W3 )L(l2,la')) e s (l 3 )6 s (l 4 



-c?* cos2(^ ll + ip h ) (L(i 2 , -i 1 )e s (i 3 )e s (i 4 )> - c? E cos2(^ ll + <p h ) (L(h, -i 2 )e s (i 3 )e s (i 4 )) 

|^C i f £ cos2( Wi - Wl )cos2( Wi + <p h ) (L(h,-h')L(h,h')e s (h)G s (h)) 



Contributions to the trispectrum from the first two terms come through the second order term in L, with the two </> 
terms coupling to s . In the last term, contributions come from the first order term of L similar to the pure lensing 
contribution to the trispectrum. 

After some straightforward simplifications, we can write the connected part of the trispectrum involving lcnsing- 
secondary coupling. We symmetrize with respect to the four indices by adding five terms corresponding to the 
additional permutations of (Zi, Z 2 , Z 3 , Z 4 ) and multiplying by a factor of 1/6. 

T est (l 1 ,l 2 ,l 3 ,l 4 ) = 

ho£(% {eg cos2(^ ll + <^i 2 )(l 3 • li)(l 4 • li) + Cf 2 cos2(^ ll + <p h )(l 3 • 1 2 )(1 4 • 1 2 ) 
+ [I3 ■ (li + I3)] [U • (la + U)] C|? 1+ i 3 | cos2(<^ ll+l3 - ip h ) cos2(^i 2+ i 4 - ip h ) 
+ [1 4 • (li + I4)] [1 3 • (1 2 + 1 3 )] C|f i+l4 | 008 2(951,+!, - (p h ) cos2(<p h+h - vi 3 ) + Perm, j 

(73) 

Note that the first two terms come from the first and second term in equation (Jt^) , while the last two terms are from 
the third term. 



(72) 
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